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^ Abstract 
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The deconfinement phase transition in large- iV c QCD is studied within the framework of an 
effective Polyakov-loop model, where the potential has a U(l) symmetry originating in the large- 
N c limit of a Z^ c -symmetric model. At the critical temperature, the shape of the effective potential 
allows the existence of Q-balls as position-dependent fluctuations of the Polyakov loop. Q-balls with 
spherical or axial symmetry are numerically obtained from the equations of motion of the effective 
model under consideration. The physical properties of these non-topological solitons (mass, charge 

H 

and size) are discussed, as well as their interpretation in terms of spinning "bubbles" , with various 
shapes, of deconfined matter surrounded by a confined environment. 
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I. INTRODUCTION 



The phenomenology related to the deconfmed phase of hadronic matter - the quarkgluon 
plasma - is nowadays a matter of intense investigations. Prom a theoretical point of view, 
the study of gauge theories at finite temperature is a challenging problem, while the QCD 
matter is or will be studied in heavy-ion collisions at RHIC, SPS, FAIR, and the LHC. 
Informations and relevant references on the topic can be found for example in pQ. 

For a given Yang-Mills theory at a nonzero temperature T, the Polyakov loop is defined 
as L(T, y) = Pe* 9 ^' drAo(T,y)^ Aq the temporal component of the Yang- Mills field and 
y the spatial coordinates. As usual, P is the path-ordering, g is the strong coupling con- 
stant and units where h = c = ks = 1 are used. The Polyakov loop is actually such that 
(L(T,y)) = (7^ 0) when the theory is in a (de)confined phase [2]. Since gauge transforma- 
tions belonging to the center of the gauge algebra only cause L(T, y) to be multiplied by an 
overall factor, it is tempting to conjecture that the confinement /deconfinement phase transi- 
tion might be linked to the spontaneous breaking of a global symmetry related to the center 
of the considered gauge algebra. In the particular case of SU(iV c ), deconfinement might thus 
be driven by the breaking of a global Z^ c symmetry [3], with the following dimensionless 
order parameter 

d> = — Tr c L. (1) 

Note that the color-averaged Polyakov loop <fi will simply be called Polyakov loop in the 
following. 

It is well known that the thermodynamic properties of pure gauge SU(3) QCD can be 
studied by resorting to an effective scalar field theory where the potential energy density is Z 3 - 
symmetric, with e.g. the form [I] U oc — a 2 |0| 2 + a 4 \4>\ 4 + a 3 (4> 3 + (ft* 3 ) , see also Refs. [SJE]- It 
has further been proposed in [7j that, for an arbitrary number of colors, the potential energy 
density should be of the Zjv c -symmetric form U = a-i |0| 2 +ct4 |0| 4 +a7v c ((f> Nc +(p* Nc ). Moreover, 
the large-iV c scaling and the temperature-dependence of the thermodynamic observables 
severely constrain the coefficients at large N c , as it has been shown in [8]. One then 
gets a potential that becomes U(l)-symmetric in the large- iV c limit. The main results of [S] 
concerning the Polyakov-loop effective potential will be summarized in Sec. [IT} 
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Polyakov-loop-inspired approaches are a precious tool to understand the QCD phase 
diagram: We refer the interested reader to e.g. Refs. [9], where the Polyakov loop is both 
coupled to quark fields and responsible for an effective potential. A direction that has been 
less studied so far is to work beyond mean field approximation and let the Polyakov loop 
fluctuate thanks to a standard kinetic term of the form d^d^ip*: Dynamical effects like 
hadronization through decay of the Polyakov loop [TQ1 [EE] or formation of plasma bubbles 
in heavy ion collisions [12J can then be studied. The plasma bubbles found in the above 
references are actually static and spherically symmetric solutions of the equation of motion 
associated to a given Polyakov-loop Lagrangian; they can be interpreted as a sphere of 
deconfined matter surrounded by a confined medium. One of the main goals of the present 
paper is to go a step further in this direction and show that spinning plasma bubbles, and 
more generally bubbles with axial symmetry, may also exist. Such bubbles can reasonably 
be imagined to form in non central heavy-ion collisions, where the incident nulcei have a 
relative angular momentum. 

It should be pointed out that there exist alternative descriptions of the pure gauge theory 
which make use of the transverse particle concept, see e.g. [T3TfT5] . Descriptions of the 
physical phenomenon based on the AdS/QCD correspondance principle are also available, 
see namely 



The paper is organized as follows : The model we use is presented in Sec. |III[ It is 
directly inspired from the large-iV c , U(l) symmetric, Polyakov loop potential presented in 
[8]. The advantage of a U(l) symmetry is that it allows in principle the existence of Q- 
balls. In general, Q-balls refer to regular, finite mass and localized classical solutions of a 
complex, self-interacting scalar field theory [T9"j |2"U] ; they form non-topological solitons of the 
underlying equations. In the present context, these solitons will come out as natural solutions 
of the Polyakov loop equation of motion. Q-balls with spherical and axial symmetry will be 



found by numerical resolution of the Polyakov-loop equations of motion in Sees. [TV] and [V] 
respectively. Their interpretation in terms of bubbles of deconfined plasma will be discussed, 
as well as their main physical properties. The various solutions obtained are compared in 



Sec. VI, while a discussion of the possible extensions of the present results to finite- N c cases 



and concluding comments are given in Sees. VII and VIII 
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II. POLYAKOV-LOOP LAGRANGIAN 



According the arguments exposed in the introduction, an effective Lagrangian for a finite- 
temperature SU(iV c ) Yang- Mills theory should involve the Polyakov loop ([I]), which is seen 
as a complex scalar field, as well as a Zjv c -symmetric potential, denoted V g . In order to 
mimic at best known results in Yang-Mills theory, V g should be such that: 

• The pressure p g = — min^(V g ) is proportional to N^T A at large N c and T in order to 
recover asymptotically the Stefan-Boltzmann limit for a free gluon gas. 

• The value of the Polyakov loop minimizing V g is iV c - independent at the dominant order 
[2l] . Numerically, |0o| = in the confined phase, > in the deconfined phase, and 
tends toward unity at very large T. 

• There exists a critical temperature T c above which the absolute minimum is nonzero, 
so that one has a first-order phase transition. T c has to be seen as a typical value for 
the deconfinement temperature in SU(iV c ) Yang-Mills theory since the deconfinement 
temperature appears to be iV c -independent up to corrections in 1/N% (22]- According 
to previously obtained results, the norm of the Polyakov loop when T — > T c + may be 
found between 0.4 and 0.6 [23l EE 



(2) 



It has been shown in [8] that the potential 
V g = N*T 4 a(T) 

with 

1{T) = 0.74-0.26 tanh 



2 _ 4 |0| 4 + ^4- [8/(T) 2 - 1](0^ + <f>* N °) 



T \ T 
2.10 I -0.60- 



T } T c 



and a(T) 



/(T) 4 V.135 ln(T/T c + 1.5) 



(3) 

has all the required features and is moreover such that it leads to numerical values of the 
pressure and of the Polyakov loop in excellent agreement with recent lattice data [231 123] , 
both for N c — 3 and in the large- N c limit. This last limit is particularly interesting for our 
purpose since, when iV c — > oo, V g becomes U(l)-symmetric and reads [S] 

VM\,T)= N?T*a(T) |0| 2 (1 - 4|0| 2 ) \<j>\ < /(T), 

+oo |0| > l(T). (4) 
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FIG. 1: Polyakov-loop potential V^(|0| c ) versus \(j)\ in the large-iV c limit (dashed line); V g is given 
by (|3j) and @. The potential used in our numerical calculations, that is ([5]), is also plotted for 
comparison with (3 = 10 and C = C\q = 27 according to ^ (solid line). 

It is worth mentioning that l(T c ) = 1/2, a(T c ) = 2/3, and that, as expected in the case of 
a first-order phase transition, the two minima are degenerate at the critical temperature: 
V g (0, T c ) = V g (l/2, T c ) = 0. The shape of the potential at T = T c is shown in Fig. [lj it is the 
typical kind of shape for which soliton-like solutions exist, as we will show in the following. 

Potential Q, becoming infinite once \<f>\ is large than some T-dependent value, is not 
very convenient for numerical applications such as those that are to be performed. Instead 
of using this last form, we will rather introduce the following potential 



V g {\<l>\,T e ) = |iV c 2 T c 4 (l0| 2 -4|0| 4 + C|0n 



(5) 



which can be seen as a more tractable approximation of Q at T = T c , with a term in \(f>\P 
introduced to mimic the "wall" of Q. Generic values of C can be considered, but the main 
emphasis will be put on the particular choice 



C = C h 



where 



4/3-16' 



(6) 



It is such that potential ^ is still zero in its two degenerate minima, namely = and 
|0| = \(p/3\- This second minimum tends toward 1/2 when (3 is arbitrarily large: Q is then 
recovered as the limit of ([5]) when (3 tends toward infinity. The potential ^ is shown in 



Fig. [I] for = 10 and C = C 10 = 27 according to @: O ne sees that, though not being 
an extremely accurate approximation of Q when > 0.5, it shares the same structure, 
and will presumably lead to qualitatively similar results, which is satisfactory for our mostly 
exploratory purpose. Note that |0io| = 0.58 instead of 1/2, but this is still acceptable if one 
considers values previously found in computations of the critical value of the Polyakov loop 

It sould be pointed out that the replacement of the potential Q by (J5| introduces an 
extra conserved quantity related to the continuous symmetry (for instance the charge Q 
defined in the following). The consequence of this will be commented in due course. 



III. THE MODEL AND THE EQUATIONS 



A. Reduced Lagrangian 

We now consider that the Polyakov loop is a dynamical complex scalar field. According 
to the suggestion of e.g. Ref. j5], we take a kinetic part of the form N^^d^d^cf)* / A, 
which has both the correct energy dimensions and the expected ^-scaling. A is the 't Hooft 
coupling. The starting Lagrangian, valid at T = T c , is thus given by 

N 2 T 2 o 

C phys = -^d^f - -iV 2 T c 4 (|0| 2 - 4|0| 4 + C|«/f), (7) 

where <p = y^ 1 being the spacetime coordinates. It is natural to further define dimen- 

sionless variables x M related to the original (physical) ones by 

"3 



— Iphys With Iphys — y, y 




so that the above Lagrangian can be replaced by 
where <fi = (f){x^) and 

v(i0D = i0i 2 -4i0i 4 +ci0r (10) 

The corresponding classical equations of motion read 

w = S = * (1 - 8 h 2 + ^r 2 ); (11) 



plus the complex conjugated equation. 



B. The ansatz for Q-balls and conserved quantities 

The standard way to construct Q-balls consists in looking for solutions of the form 

= exp(iux°)^(x), (12) 

where ip(x) is a function of the space variables. The real parameter u constitutes an essential 
characterisation of the stationary solution. As will be soon discussed, solitons only exist for 
ui taking values in a finite interval. 

Two types of solutions will be searched for: Spherically- or axially-symmetric ones. 
Spherically symmetric solutions correspond to further assuming i]){x) = x( r )> where r is 
the radius in spherical coordinates. Axially symmetric solutions correspond to the ansatz 
if){x) = exp(ikip)x(r, 9), where {r, (p, 9} denote the standard spherical coordinates and where 
k is an integer. We remark that x is assumed to be a real function. 

The solutions we will build can be characterised by their energy M and by a dimensionless 
conserved charge Q, respectively defined by 

M = M phys J d 3 x Too (13) 

and 

Q = 2co [ d 3 x |0| 2 . (14) 



The temporal component of the energy-momentum tensor represents the energy density, 
given by 

T oo = w 2 |0| 2 + V0-V0* + \/(|0|). (15) 
The natural mass scale introduced above reads 

2A^" (16) 
The conserved charge Q finds its origin in the U(l)-symmetry of the considered lagrangian, 
leading to a conserved Noether current of the form J M = i^d^cj)* — 0*6^0), Q being the the 
space integral of J . Note that the axially symmetric solutions having k ^ are spinning 



Q-balls whose angular momentum J is related to the charge Q according to J = kQ [26J. 
This relation is specific to Q-balls: the solutions constructed with the potential (J2j) could 
still have a conserved angular momentum while Q would be meaningless. 

Let us point out that all the masses and lengths to be plotted below are expressed in 
units of M phys and l v h ys respectively. 

C. Physical quantities 

At this point we have to stress that the proposed model does not intend to provide an 
accurate description of "realistic" plasma bubbles that might form in heavy ion collisions. 
It is rather a first step to explore the different kinds of soliton-like solutions of a typical 
Polyakov-loop-inspired model of the QCD phase transition. Nevertheless, it has been shown 
in Ref. [8] that our model correctly reproduces the QCD equation of state in the large- N c 
limit computed on the lattice [25J, and leads to a commonly accepted shape for the QCD 
phase diagram once coupled to quarks. Therefore we think that the Q-balls we will obtain 
have something to tell, at least qualitatively, on non-trivial phenomena that might arise at 
the deconfinement temperature. 

It is worth estimating here the physical mass and length used in the model thanks to some 
known numbers concerning the quark-gluon plasma at N c = 3. First, a typical value for the 
deconfinement temperature in QCD is T c = 0.2 GeV []]. Second, the 't Hooft coupling can 
be expressed as a function of the strong coupling constant a s through A = 4N c 7ta s . A way 
to estimate a s is to note that the short-range part of the static interaction between a quark 
and an antiquark scales as — (4/3)a s /r , at least from T = to T c . Recent lattice studies, 
performed at N c = 3 and with 2 light quark flavors, favor a s = 0.2 up to T = T c [34] , that 
is the value we retain here. We are then in position to estimate that, at N c = 3, 

l phys = 2.2 GeV" 1 = 0.44 fm, M phys = 0.107 GeV. (17) 
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IV. SPHERICALLY SYMMETRIC SOLUTIONS 



Starting from Lagrangian (10), the field equation reduces, in the spherically symmetric 
case, to the differential equation : 

which has to be solved for r e [0,oo]. The regularity of the solution at the origin implies 
dx{0)/dr = 0, the fmiteness of the energy and the charge impose x(°°) = 0. These conditions 
specify the conditions at the boundary 

Adapting the results of [26J to our potential, it can be shown that Q-balls with spherical 
symmetry exist for 

2 



4/3 — 16 / 8 \ 
u 2 p < u 2 < 1 with uj 2 p = min{0, 1 - -^—^ \^ —-^ \ } (19) 

In particular, up = if C < Cp. Only for C > Cp the lower bound is non zero. The case 
mostly studied in [26j 12Z] corresponds (up to an appropriate renormalisation of the field and 
of the radial variable) to /3 = 6 and C = 44/10 so solutions exist for 1/11 < oo 2 < 1. In 
the present work, we are mostly interested in the case C = Cp for which Q-balls exist for 
cu G [0, 1]. For = 6, i.e. the most commonly used value in studies devoted to Q-balls, we 
have Cq = 4. 

In order to calibrate our solutions with cases studied in the literature, we have first solved 
the equations for C = 44/10 and several values of (3. The mass M and charge Q of these 
solutions are reported as functions of to on Fig. [2] (left side). On the right side, the same 
quantities are reported for N = 10 and two values of C : the value C = 44/10 corresponds to 
Ref. [26J, and C = 27 = C\q is the case under consideration for this paper. We have checked 
that we recover indeed previously obtained solutions for the aforementioned values of the 
parameters. 



It is worth saying that Eq. ( 18 ) has, to our knowledge, no analytic solution for arbitrary 
values of /3 and has to be solved numerically. We exclude of our considerations the trivial 
case C = f3 = 4, actually leading to a free Schrodinger equation for \- It should be pointed 
out that even a numerical approach is not obvious because the boundary conditions are 
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FIG. 2: Left: Mass (solid lines) and charge (dashed lines) of spherically symmetric Q-balls for 
C = 44/10 and (from top to bottom) /3 = 6 ,8, 10. Right: Idem for f3 = 10 and (from top to 
bottom) C = 27, 44/10. 

compatible with the trivial solution xi r ) = 0. As a consequence, obtaining a non trivial 
numerical solution requires a starting guess which is reasonably close to the desired solution. 
We used a collocation method for boundary-value ordinary differential equations, equipped 
with an adaptive mesh selection procedure |28j. Our solution were constructed with a relative 
error of order 10~ 8 . 

In the following, we focus on the case = 10 and C = C\q = 27. From a technical 
point of view, the numerical construction of the solutions, especially the axially symmetric 
ones, becomes more involved while increasing 0; that is why no values of large than 
10 will be used. Nevertheless we think that the case /3 = 10 already captures most of 
the qualitative features of the limit /3 — > oo, that corresponds to our original Lagrangian. 
Solutions corresponding to w = 0.1, 0.3, 0.6, 0.95 are plotted in Fig. [3] in order to show the 
behaviour of xi r ) an d T 00 (r) for various values of u lying in the range [0, 1]. 

For small values of u, the minimum of the effective potential is deep and the scalar 
field "likes" to stay close to this minimum: It can be observed in Fig. [3] that, confirming 
this observation, our numerical results demonstrate that the radial function x(r) remains 
practically constant (i.e. x( r ) ~ x(0) = I'Ppl) inside a large sphere centred at the origin. 
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FIG. 3: Profile of x( r ) (solid lines) and of the energy density (dashed lines) in the spherically 
symmetric case for (3 = 10, C = C\q = 27 and (from left to right) u = 0.95, 0.6, 0.3, 0.05. 

Then it brutally decreases to the the asymptotic value x{ r ~~ °°) = 0, forming a thin wall. 
Accordingly the energy density presents a plateau for T o = w 2 0/3 inside the sphere, then 
reaches a maximum on a spherical shell corresponding to the wall and finally decreases to 
zero. These features appear on Fig. |3j The scenario is quite different for u> close to unity. 
Here the effective potential's minimum is not deep anymore and the scalar field quickly 
reaches its asymptotic value as shown in Fig. [3j This is the thick- wall limit. 

A question that is worth asking is: Does the Q-balls we have found describe some de- 
confined sphere - or spherical shell - surrounded by a confined environment? Although the 
Polyakov loop is a straightforward indicator of confinement in a mean-field treatment, this 
is less obvious in the present framework, where it is never exactly zero, unless asymptoti- 
cally by construction. We find that a good criterion could be instead given by the value of 
the energy density. According to the recent lattice computations of the pure gauge QCD 
equation of state given in [25], one can estimate that the deconfined phase is characterised 
by energy densities such that e > 0.35N%T£. By recalling that e = To M p hys/lp hys , we can 
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translate the criterion e > 0.35 into 

Too > 0.45. (20) 

This last inequality can be used in the following as a criterion to seprate the confined regions 
from the deconfined ones in the solutions we find. A glance at Fig. [3] shows that only Q-balls 
with uj > 0.6 show deconfined regions according to the above criterion. Only these Q-balls 
could then be interpreted as "plasma balls" or "plasma shells" , while the others would just 
be fluctuations of the Polyakov loop in a completely confined region of space. As it can be 
observed in the left panel of Fig. [2j Q-balls with u around 0.9 not only contain a deconfined 
part, but are also the lightest ones that can be obtained. The Q-ball with minimal mass 
is reached for u = 0.88, with a charge Q = 13.7 and a mass M = 14.7 (1.57 GeV). In the 
thick-wall limit, that is for u = 1, one has finally Q = 55.0 and M = 55.2 (5.87 GeV). We 
mention finally that the Q-ball with oj = 0.6 has a charge Q = 44.5 and a mass M = 35.5 
(3.78 GeV). 



V. AXIALLY SYMMETRIC SOLUTIONS 



Let us now turn to axially symmetric solutions. To our knowledge, they have never been 
obtained so far within the framework of a Polyakov-loop model or, more technically, in a 
Q-ball study using the value (3 = 10. The field equation reads in this case 

d 2 X 2 dx 1 d 2 x cos6 d X k 2 = _ 2 R 2 PC . 2) 

dr 2 rdr r 2 d9 2 r 2 sine d9 r 2 sm 2 e X X{ U X 2 X h 1 ] 

and has to be solved in principle for r G [0, oo] and 6 G [0, 7r]. The above equation is noto- 
riously difficult to solve because it is an elliptic non-linear partial differential equation. It 
has to be completed by appropriate boundary conditions which are given below. However, 
rendering the numerical construction more difficult, the different sets of conditions are com- 
patible with the trivial solution x(r, 6) = 0. Several solutions of equations of the type above 
were constructed in different contexts, see e.g. [29TI3T] . Comparisons of these solutions with 
ours constitutes a useful crosscheck of our numerical method. For our problem, the integra- 
tion was performed numerically by a routine [32] based on the Newton-Raphson method. 
Concerning the angular dependence of the solutions, we will focus on Q-balls which are even 
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or odd under the reflection — > n — 0. Accordingly the interval of the variable can be 
limited to G [0,7r/2]. In addition to their mass and charge, the solutions with k > are 
further characterized by an angular momentum given by J = kQ (see [26] for details). 



A. Spinning solutions 



Setting k — 1, Eq. (21) can be solved by imposing the boundary conditions 



X (O,0) = O , X (oo,0) = O , x (r,0 = O) = O , d eX {r,0 = tt/2) = 0. (22) 

Solutions of this type are even under the reflexion n — and have angular momentum 
J = Q. Such Q-balls have been first obtained in |26j for /3=6; they can be deformed to the 
case of our potential. As an illustration, plots of the function x(r, 0) and of the region where 



the energy density satisfies to the criterion (20) are presented in Fig. 4 for the solution 



u> = 0.8 . It is readily seen that the spinning Q-ball found by applying the boundary 



conditions (22) has a toric shape, winding and spinning around the z-axis. Such a Q-ball 
may be seen as a "plasma ring" : A spinning torus of deconfined matter. Note that all the 
values of u for which we have found a solution show a similar toric shape and a deconfined 
part. 

A more accurate description of the \ field and of the energy density in the equator plane 
is presented in Figs. [5] for several values of u. We were able to construct solutions for 
u G [0.55,0.99]. The upper limit for u is the same as in the spherically symmetric case as 
it has previously been observed in [27]. It is likely that solutions exist also for smaller u 
but their construction become highly involved and understanding in detail the full pattern 
of solutions is not the aim of this exploratory work. The numerical difficulties in obtaining 
solutions for small u are recurrent and also pointed out in [27]. In the case of spherically 
symmetric solutions, it turns out that several solutions exist with the same 0(0) and different 
u, the values of 0(0) becoming weakly dependant on u. In the case of axial solutions, the 
difficulty seems to be related to the fact that the solution spreads over spacetime and the 
mesh has to be updated constantly. At least, the spinning solution with the minimal energy 
occurs inside the interval that we explored and is reached for u> m 0.86 where M ~ 55 (5.9 
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FIG. 4: Left: Plot of the scalar field function in the = plane for the spinning Q-balls with 



boundary conditions (22) corresponding to @ = 10, C = C\q = 27 and u = 0.8. Right: Plot of 
the region where Too > 0.45 for the same solution. Note that {x, y, z} denote standard cartesian 
coordinates. 

GeV) and Q fa 51. The largest diameter of the deconfined part in this lightest solution is 
about 8, that is around 3.5 fm. 



B. Further axial solution 



The angular dependence of the spherically symmetric and spinning solutions discussed 
above have the same symmetries as the y o ° and oc sin(0) spherical harmonics. It is 
therefore tempting to assume that families of solutions exist with the symmetries of the 
Y£(9,<p) spherical harmonics with — m < k < m. For instance the parity-odd spinning 
solution of |26j has the same angular symmetries as the function Y%. According to this 
observation there should exist also within our framework a family of solutions related to 
e.g. the function Y® oc cos(#). Such solutions were first constructed in [33] with (3 = 6 
and different values for the parameter C, and here checked to exist for generic values of the 
potential's parameters. 
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FIG. 5: Left: Energy density in the equator plane for the spinning Q-balls with boundary conditions 
( |22[ ), N = 10 and C = C\q = 27 for several values of u. From top to bottom: u> = 0.8, 0.7, 0.6 
(solid lines) and w = 0.9, 0.95, 0.99 (dashed lines). Right: Idem for the field x( r ) 7r /2) with, from 
left to right, u = 0.8, 0.7, 0.6 (solid lines) and u = 0.9, 0.95, 0.99 (dashed lines). 

Configurations with the Y® angular dependence can be obtained with j3 = 10 and C = 27 



as Q-ball solutions of Eq. (21) and the boundary conditions 



X(O,0) = O , x(oo,#) = , d X (r,9 = 0) = , X (r, = tt/2) = 0. (23) 

These boundary conditions ensure the solution to be odd under parity. We are able to 
construct with a good accuracy the branch of solutions for ui G [0.55,0.95]. On this interval 
a configuration with minimal mass seems to be reached for ui « 0.92, corresponding to 
M ~ 26 (2.8 GeV). As an illustration, the solution obtained for uj = 0.8 is plotted in Fig. 
[6| The norm of the x -ne ld and the energy density are concentrated in two regions of the 
z-axis, say around z = ±z c : We find z c w 6 (2.6 fm) for u = 0.8. Such Q-balls could be seen 
as a two-center deconfined region, much like a bound state of two "plasma bells" . 

The evolution of the energy density on the (positive) z-axis for different values of u can be 
estimated from Fig. [7j It suggests that the structure of the lump may become quite involved 
when to decreases. Again, the construction of such low-w Q-balls has not been undertaken, 
and deconfined regions appear for all the considered values of uj. 
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FIG. 6: Left: Plot of the scalar field function in the = plane for the Q-ball with boundary 



conditions (23) corresponding to j3 = 10, C = Cxo = 27 and oj = 0.8. Right: Plot of the region 



where Too > 0.45 for the same solution. Note that {x, y, z] denote standard cartesian coordinates. 

VI. COMPARISON AND STABILITY 

To complete the description of the different solutions just found, we report their mass 
and charge as functions of the parameter u on the left panel of Fig. 8. This part of the 
figure is limited to u G [0.6,0.95], where the three types of solutions have been constructed. 
Likely due to the large exponent of the field /3 — 10, the numerical resolution of the partial 



differential equation (21 ) becomes unreliable for uo < 0.6; a more appropriate discretisation of 
the space of integration would probably do the job. However, as stated above, we have been 
able to obtain the solutions with the lowest energy, that are the most likely to be produced 
at the deconfinement phase transition. Let us label the different branches of solutions by 
the couple (m, k) according to their angular dependence in relation with the Y£ spherical 
harmonic. As expected, for fixed uj, the solution with the lowest energy is the spherically 
symmetric solution, then come successively the axially symmetric (1,0) and the spinning 
(1, 1) solution. 
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FIG. 7: Field (solid lines) and energy density (dashed lines) on the positive z-axis for the Q-balls 
with boundary conditions (23), (3 = 10 and C = C\q 



27 for different values of oj: From left to 
right, oj = 0.6, 0.75, 0.9, 0.95. The corresponding values on the negative z-axis may be obtained 
by recalling that x (Too) is an odd (even) function of z. 

Perhaps more relevant is the right panel of Fig. 8 which represents a (Q, M) plot of the 
different solutions. As observed in [27] for example, it reveals that, for each family (m,k), 
the curve M(Q) presents different branches terminating into spikes at critical values of the 
charge. According to standard arguments based on bifurcation theory, the branch with the 
lowest energy should be stable while the branches with higher energy are unstable, see e.g. 
[35]. 

For all values of the conserved charge covered by our solutions, the solution with the 
lowest energy is a spherically symmetric solution. On the graph, this corresponds to the 
lowest branch. Note that the two lines corresponding to the spherical solution are plotted 
up to Q = 200 on the graph, but actually terminates with a spike at Q ~ 2200, M m 806 
(the graph was cut for obvious reason). 

For the values Q > 55, the two axially symmetric solutions coexist and it can be checked 
that the spinning solution always possesses the lowest mass at a given charge. In the region 
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FIG. 8: Left: Mass (solid lines) and charge (dashed lines) as function of w of, from top to bottom, 
the spinning solution (1,1), the axial solution (1,0), and the spherical solution (0,0). Right: Mass 
versus charge of the different solutions obtained. A zoom on the region with the lowest masses and 
charges has been also added. 

27 < Q < 55, the spinning solution does not exist and the second lowest mass is the (1,0) 
axial solution (corresponding to the red line). More details of this branch are shown in the 
window, which is a zoom of the plot on the low Q and M region. Comforting these results 
about the stability, obtained on the basis of the theory of bifurcations, let us finally point out 
that the various branches of stable solutions fulfil the condition M < Mf ree , Mf ree = m^Q 
where is the mass of the boson in the underlying field theory (with our units = 1) 
and Q is the charge which is, in conventional Q-ball literature, is interpreted as the particle 
number. The Q-balls located under the M = Q line in the right panel of Fig. 8 are then 
stable. 

Several families of solutions could in principle be constructed extending the above results. 
Let us mention three possibilities : (i) radial excitations of the spherically symmetric solution, 
where the scalar field presents nodes at some values of the radial variable (see e.g. [261 [27]): 
(ii) further angular excitations corresponding to higher values of the integers k and m; (iii) 
a mix of (i) and (ii). 
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VII. FINITE N c 



The large-iV c limit leads to an effective potential with continuous U(l)-symmetry and 
consequently to Q-balls. Although a detailed numerical investigation of this case is out the 
scope of the present paper, it is worth commenting the situation at finite N c . 

When N c is finite (and > 2), the symmetry of the Lagrangian is a discrete, Zjv c one 
enforced by the potential term aN c (4> Nc +4>* Nc ), causing the appearance of a term aM c N c (j) Nc ~ 1 
(a^iVc^* 7 ^ -1 ) in the equation of motion for <p (4>*). One readily sees that the ansatz (12) 



leads to incompatible equations unless u — 0. In this case, the equations of motion (18) and 



(21) are recovered. So all the Q-ball solutions with oj = existing at infinite N c are expected 
to exist at finite N c too. In particular, we have found such solutions in the spherically 
symmetric case, and we have no argument suggesting that axially symmetric solutions with 
vanishing u do not exist, although they are technically complicated to build. 

Which quantum numbers could be used to label these finite- N c solutions ? The charge Q 
labelling Q-balls is a consequence of the U(l) symmetry, so it is no more relevant. However, 
Q is linked to the total angular momentum J = kQ of the solution, which is defined from 
the energy-momentum tensor, as well as the mass M of the solution. So M and J are still 
good observables to classify solutions at finite N c , instead of M and Q. 

The Q-ball-like configurations that one could find at finite N c are likely to be unstable 
because of the absence of a conserved charge Q. From the paper [TO], we already have 
evidences that N c = 3 configurations of the form 4>(x°, x) = a(x°) exist, where a(x°) start 
from the "deconfined minimum" of the potential and falls down to the confined one as 
time increases. This is an example of a solution modelling an unstable configuration of 
deconfined matter. So we think that unstable solutions of the type 4>(x°,x) may exist, such 
that 0(0, x) = x{%) can be identified with of the aforementioned uj = solutions, but with 
4>(oo,x) = 0. This kind of solutions could be interpreted as an unstable plasma ball. 
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VIII. CONCLUSIONS AND OUTLOOK 



In this paper, we have used some knowledge of Q-balls developed several years ago in a 
context of non-topological solitons within an effective model attempting to describe quali- 
tatively the deconfming phase transition in the large- N c limit. The effective theory used is 
a scalar field theory with U(l)-symmetric potential (seen as the infinite N c limit of a Zn c ~ 
symmetric theory), where the complex scalar field is the color-averaged Polyakov loop. At 
the deconfinement temperature, the effective potential has two degenerate minima mimick- 
ing the first-order deconfinement phase transition. Such an effective potential has a shape 
that strongly resembles to potentials typically used in the study of Q-balls, motivating the 
present study. 

For technical convenience, the infinite wall of the effective Polyakov loop potential oc- 
curring for large norms of the effective field have been replaced here by a large power of 
the scalar field. As expected, the main families of Q-balls, obtained in previous studies, 
can be constructed for this potential too. Through an investigation of the energy density 
distribution in the different Q-balls obtained, we have been able to interpret some of them 
as bubbles of deconfined matter, or quark-gluon plasma. The different shapes we have found 
correspond to (by increasing mass for a given charge): spheres, tori, or two-bell bound states. 
These last two solutions have been obtained here for the first time within a Polyakov-loop 
model. The physical masses and sizes of the obtained solutions cover a large range of values; 
we recall that the lightest Q-ball is a spherically symmetric one with a mass M « 1.6 GeV 
and and a typical size of 1.5 fm for the deconfined region. Torus- and two-bell like solutions 
are systematically heavier (and larger) than this solution. 

Returning to potential ([2]), it could be mentioned that our Q-balls corresponding to to = 
could provide solutions of the more realistic, finite-A c , Zn c equation. We believe that our 
Q-balls with generic values of u could produce mean values (or suitable initial profile) of 
more realistic time-dependent solutions of the potential (|2]). As far as the stability of the 
solutions is concerned, the non conservation of Q for the realistic equations would allow the 
lump to decay, as expected in the context of the fireball expansion in hadronic collisions. 

The precise behaviour of the solutions in regions of low u were not pushed to details since 
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it was not the aim of this exploratory study: A more detailed classification of solutions, with 
higher angular momentum in particular, as well as a more realistic description of the gluon 
plasma - through the building of unstable, finite- N c , solutions for example- is left for future 
works. 
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